BGU-PH-95/06, CERN-TH/95-144 



Relic Gravitational Waves from String Cosmology 

R. Brustein( a,c ), M. Gasperini^, M. Giovannini( b ' c ) and G. Veneziano^ 

^Department of Physics, Ben-Gurion University, Beer-Sheva 84105, Israel 
( b > Dipartimento di Fisica Teorica, Via P.Giuria 1, 10125 Turin, Italy 
( c ) Theory Division, CERN, CH-1211, Geneva 23, Switzerland 



Abstract 

A large class of string-cosmology backgrounds leads to a spectrum of relic stochastic grav- 
itational waves, strongly tilted towards high frequencies, and characterized by two basic 
parameters of the cosmological model. We estimate the required sensitivity for detection of 
the predicted gravitational radiation and show that a region of our parameter space is within 
reach for some of the planned gravitational- wave detectors. 
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1 Introduction 



It is notoriously difficult to find accessible experimental signatures of fundamental strings 
because of their small Planckian size. Possibly, an interesting exception to this rule is rep- 
resented by the cosmological predictions of string theory, since these originate from physics 
of the Early-Universe, when space-time curvatures may have been of Planckian strength. 

In order to arrive at some concrete, yet generic, predictions of string cosmology, we shall 
consider a large class of models in which a period of dilaton-driven inflation |], @, |3[] is 
followed by a stringy epoch, during which the curvature remains of the order of the string 
scale X~ 2 , and then finally, after possibly a short dilaton-relaxation era, by the standard 
(radiation then matter dominated) cosmology. As discussed in detail in the presence 
of a high-curvature stringy epoch appears to be unavoidable for a viable inflationary string 
cosmology scenario. Additional support to this point of view was given in ||. 

Recently, in collaboration with V. Mukhanov M, we have discussed the main properties 
of metric perturbations in a dilaton-driven background. Particular attention was given to 
the correct treatment of scalar perturbations which, in a standard treatment, appear to grow 
too large for the applicability of linear perturbation theory. By carefully "gauging down" 
certain growing modes, we were able to show that both scalar and tensor perturbations can 
be treated perturbatively and that they exhibit very similar spectra. A characteristic feature 
of these spectra is that, unlike the spectra of the standard inflationary scenarios, they are 
not flat but strongly tilted towards higher frequencies, as originally noted in [[/]]. 

In this note we will concentrate on tensor perturbations, those associated with gravita- 
tional waves (GW), at frequencies which may be accessible to earth-based GW detectors. 
Quite amazingly, their spectrum turns out to be rather independent of most of the details of 
the above string cosmology scenario. As shown in the first part of this paper, the spectrum 
can be completely given in terms of two parameters, the value g s of the string coupling pa- 
rameter at the end of the dilaton-driven phase (equivalently, at the beginning of the string 
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epoch), and the total red-shift z s occurring during the string era. The theoretical advantage 
of considering GW signals stems from the fact that, unlike the electromagnetic perturba- 
tions which underwent a complicated history until recombination, gravitons decoupled since 
right after the Planck era. As a consequence, their present spectrum should be a faithful 
portrait of the very early Universe. On the other hand, the detection of a cosmological 
background of GW requires extremely precise length measurements, typically at least of the 
order of 5L/L < 10~ 21 ||. In the second part of this paper we will show that the expected 
range of g s and z s includes a region which should be accessible to future gravitational wave 
experiments. 

Throughout this paper we shall be working in the so-called String-frame, in which weakly 
interacting strings move along geodesic surfaces. Identical results would follow by adopting 
the more conventional Einstein frame in which the curvature is canonically normalized, but 
we believe the physics to be more transparent in the former frame. In the String frame the 
string length parameter which is the short-distance cut-off of string theory, X s = \/ a'h, is 
constant, while the Planck length, Ap = y/Gjyh, evolves in time as Ap = e^/ 2 A s in a time- 
dependent dilaton background. In our scenario, the background evolution starts from the 
string perturbative vacuum ||, therefore Ap grows from a very small initial value to a value 
g s X s reached at the beginning of the stringy era and, finally, to its present (very large!) value 
of about 10~ 33 cm at the beginning of the radiation dominated era. 

2 Primordial gravitational- wave spectra 

Let us consider, for the moment, an isotropic, (3 + l)-dimensional, spatially flat cosmology. 
Following (see also 0) it is easy to show that the canonical variable ip associated with 
tensor perturbations is related to the String-frame metric g^ v via 

where a(t) is the isotropic scale factor, rj^ is the flat Minkowski metric and g = exp(0/2). 
The Fourier modes of each of the two physical, transverse-traceless, polarizations satisfy the 
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following simple wave equation 



V£ + [k 2 - V( V )]ij k = 0, 



V( V ) = (g/a)(a/g)" 



(2.2) 



where a prime denotes differentiation with respect to conformal time rj (adrj = dt) and k is 
the comoving wave number, related to the physical one, u, by k — ua. Note that, since a/g 
has power-like behavior in i] during the dilaton-driven phase |l], @, §], ^(r/) grows like i]~ 2 
during that epoch reaching a maximum at 77 = r] s , i.e. when iJ^ 1 ~ 7] s a s ~ A s = M" 1 . We 
expect ^(77) to keep growing during the stringy era and then to fall rapidly to zero at the 
onset (77 = rji) of the radiation-dominated era since, in that phase, a/g ~ 77. 

A given mode k will be well inside the horizon initially, then hit the potential barrier 
V(r)) at some "exit" time r] ex ~ k" 1 , and leave the barrier at some later "reentry" time 
V = Vre ~ Vi- The approximate solutions of eq. (|2.2j ) in these three regimes are given by 

(2.3) 

(2.4) 
(2.5) 



A k 



T] ex <T] < T] re 



Vk 



c + (k)e~ ikri + c-{k)e ikr) 



V > Vr 



Equation (|2.3|) enforces the proper normalization of the primordial vacuum fluctuations. In 
the regime described by eq. (|2.4|) the perturbation is frozen outside the horizon and the 
two terms appearing there correspond to the freezing of h and of its canonically conjugate 
momentum, respectively. Finally, in Q2.5p , the magnitude of c_ (so-called Bogoliubov coeffi- 
cient) gives the amplification of the GW with respect to a minimal vacuum fluctuation. The 
actual value of |c_| can be easily obtained by matching the above solutions and their first 
derivatives at each transition time, 

2 



2kr/ 1 \c-(k)\ 



9ex/ ®e 

g re j (x r 



1 + krj! ( 9re ^ re \ + k(g ex /a ex ) 2 j drj {g/af 



(2.6) 



At this point we have to insert some information about the background evolution. In the 
simple case at hand of a D = 3 + 1 isotropic cosmology with static extra dimensions, the 
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dilaton-driven inflationary background is simply given by [[| 

a(r]) = (-r/r 1 ^, (j){rj) = -V3\n(-7]) , a/g ~ (-r?) 1/2 , -oo < rj < (2.7) 

while, for the string era, we will assume that H and d t (j) are approximately constant and 
of order A7 1 . Limiting our attention for the time being to those scales which crossed the 
horizon during the dilaton-driven phase, we thus arrive at the following estimate 

2k Vl \c4k)\ ~ (k/k s ) 1/2 z s (g s /g 1 ) [l + zf{ gi /g s ) 2 + \n(k s /k) + i] , k < k s (2.8) 

where we have denoted for convenience a re /a s ~ a\/a s = z s . In eq . ( |2~8|) k^ 1 ~ r/ s ~ (ff s a s ) _1 
is the last scale exiting during the dilaton-driven phase, gi = exp(0i/2) is a number of order 
unity which may be determined in term of the (known) present value of the ratio X p / \ s [|T(J , 
and / is the fc-independent quantity 



w/vs drj (g/a 



,2 



l + zf{ gi /g s f (2.9) 



Vs (g s /a s ) 2 

The r.m.s. perturbation amplitude over a comoving length scale k~ l (see for instance U) is 
given, in general, by \Shk(f]) \ ~ fc 3//2 |/ifc| = (g /a)k 3 ^ 2 \-ipk\- For i] > r q 1 , we then find 



\xu 1 \ 1 Hiai ( k \ 9s 
\8h k (r})\ ~ — — z s 

a{r))M p \k s J g x 



2 



(2.10) 



Equation ( p,10|) is the main result of this section and we will use it subsequently to estimate 



the required sensitivity for detection. 

It is worth stressing that the leading term in |c_| comes from the integral appearing 
in eq. ( |2.4|) , associated with the freezing of the momentum variable conjugate to h. This 
unusual result is due to the fact that, in the Einstein frame, the scale factor is = a/g, 
and our background corresponds to a contracting, rather than expanding, Universe. The 
amplification of tensor perturbations in a contracting Universe was first considered long ago 
in 
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In spite of the presence of a high-curvature regime, we expect our estimates to be valid 
for scales that went out of the horizon in the dilaton-driven phase, since they follow from 
the general physical principle that a perturbation and its canonically conjugate momentum 
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should remain frozen while outside the horizon. The perturbations thus evolve in a purely 
kinematical way, giving rise to the logarithmic term in ( |2.8| ) from evolution during the 
dilaton-dominated phase and to the second term / from the stringy epoch. 

Finally, let us digress a moment to show the stability of the result (|2.10|) with respect 
to 0(d,d) transformations which connect |13[] different homogeneous string cosmologies. A 
simple derivation of this nice feature is obtained by working in cosmic time and by using 
directly the amplitude h [defined in eq. ( j2.1| )] rather than the canonically normalized per- 
turbation if). It is straightforward to check that, in arbitrary 0(d, d) -related Bianchi I- type 
backgrounds (including possible dynamical internal dimensions and an antisymmetric tensor 
B^v), the Fourier modes of h satisfy the following simple equation (see also 0) 

K - tK, + u 2 K = o (2.11) 

where dots denote derivatives with respect to cosmic time and = <p — In \ detg^ v \ 1 ^ 2 is the 
so-called shifted dilaton, invariant under 0(d,d) transformations. Using the fact that, in 
any dilaton-driven vacuum cosmology, ~ — lnt, we easily obtain, for perturbations well 
outside the horizon, 



h u (t) ~ In 



t 


~ In \ujt\ 


tex 





(2.12) 



showing that the spectrum of GW is independent of the chosen string-cosmology background, 
and increasing our confidence that eq. (|2.10| ) is indeed the generic GW spectrum of a 
large class of string cosmology models. Note, incidentally, that this is not the case for the 
electromagnetic perturbations discussed in [p 



3 Observability 



In order to discuss the observability of our signal it is useful to rewrite our main result ( 2.10|) 
in terms of present, red-shifted proper frequencies u = k/a. One easily finds 

2" 



\5K 



M -Zeq 9sZ s 



u s = k s /a ~ z-^yiloMsZ; 1 = z^uj x ~ zj 1 g 1 ^ 2 10 11 Hz 



1 l/2 1n H; 



(3.1) 
(3.2) 



where z eq = a/ a eq ~ 10 4 takes into account the transition from radiation to matter dominance 
at t — t eq , u>i = Hidi/a ~ 10 11 Hz is the maximal frequency reached during the string phase, 
M s = Xj 1 ~ Hi, and Hq ~ 10~ 18 Hz is the present value of the Hubble scale. 

It is also convenient to rewrite our results in terms of another commonly used quantity, 
the fraction of critical density, flow = Pgw/Pc stored in our GW per logarithmic interval 
of uj. Defining dQcw/id^uj) = uj a \c-\ 2 / (M P H) 2 we have 



GW 



Z eq 9s I . . 



,9s 



uj 



jr) \SK 

JtLq, 



(3.3) 



din a; 

It emerges from eq. ( |3. 3|) that u s plays the role of an effective temperature in the sense 
that, below uj s , the spectrum is Planckian (up to logarithms of uj). The normalization of the 
spectrum, however, is different from Planck's because of the further amplification due to the 
stringy phase. Also, we do not expect the spectrum to stay Planckian above uj s (see below), 
but rather to keep growing and to reach its maximum at ui before falling exponentially. 

Finally we give, with the appropriate caveats, the generalization of the above results to 
frequencies whose exit occurred during the stringy phase, i.e. to the high frequency part 
u s < uj < uji. We are aware of the possible dangers in using field theoretic methods to discuss 
perturbations in this regime. However, in the absence of a full string theoretic calculation, 
we shall present our results as an indication of what a possible outcome might be. One finds, 
after straightforward calculations, 

uj \ 2 -P 



\Sh„ 



GW 



din a; 



0n 



a 2 z~ l 



'-HQ -1/4 



UJl 



-) 

UJl J 



6-2/3 



+ 



UJl J 



2ff 



UJl 



UJ S < UJ < UJl 



(3.4) 



(3.5) 



where j3 = — \og{g s / g\)J log z s is also the average value of g/(gH), which we have assumed 
to vary little during the string phase. We have also used the fact that during the string phase 
the curvature stays controlled by the string scale A s so that, in the String frame, the metric 
describes a de Sitter-like expansion with z s = ai/a s = r] s /r]i (see [[HJ for further details, and 
for a different derivation of the same spectrum in the Einstein frame). 

We would like to discuss now the prospects of observing our spectrum in gravitational 



wave detectors. Our main emphasis will be on the planned large interferometers LIGO Wfi 
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and VIRGO |L7|, which are expected || to start operating at sensitivities (for detection of 
a stochastic GW background) of dQcw /d In u = 10~ 6 in a frequency band around a few 
hundred Hz, and have set the ambitious final sensitivity goals of dflcw /d In u = 10 -10 in 
a frequency band around uj = 100Hz. It may well be, especially in the first stages of 



operation, that coincidence experiments between bars and interferometers |L8] could also 
be able to reach similar sensitivities at frequencies around 1 KHz. We will mention later 
other possible devices which seem to have some good potential sensitivity, especially in the 
higher frequency range, but which have not yet matured into concrete operating systems. 
Detection of stochastic GW backgrounds at frequencies below 1 Hz does not seem accessible 
with current technologies, and we therefore limit our attention to the range above 1 Hz. 

We are interested in finding the regions in our {z s ,g s } parameter space that may be 
accessible to experimental detection. From eq. ( |3.2| ) we can immediately see that the 
accessible region requires large values of z s . We may distinguish values of z s in the range 
z s < 10 9 (i. e. uj s > ujj), in which the observable spectrum at uj comes mainly from 
perturbations that crossed the horizon during the dilaton-driven era, from those in the range 
z s > 10 9 (u s < ujj) in which the observable spectrum comes mainly from those perturbations 
that crossed the horizon during the stringy era. The predictions in the range z s > 10 9 
should be considered as less robust than those in the range z s < 10 9 . In addition, we may 
distinguish values of g s in the range g s ^gi in which the dilaton does not change much during 
the stringy phase, from those in the range g s <C g\ in which the dilaton changes by a large 
amount during the stringy epoch. 



In all cases we have to impose the bound following from pulsar-timing measurements [19 



which implies dflcw / 'd\nuj<10 6 at up = 10 s Hz. We also accept the bound flcw^O.l, 



imposed by standard nucleosynthesis pOfl . Moreover, we have derived the spectrum in the 
linear approximation, expanding around a homogeneous background. We have thus to im- 
pose, for consistency, that the amplified perturbations have a negligible back-reaction on 
the metric, namely dQaw /d In uj < 1 at all frequencies and times. This, together with the 
nucleosynthesis bound, requires a range of parameters corresponding to a spectrum which is 
growing also in the stringy phase, < (3 < 3, namely (g s /gi) < 1 and (g s /gi) > z s 
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Inserting the appropriate numbers in eqs. ( |3.3|) , ( |3.5|) we find, respectively, the following 
conditions for detectability in interferometers, i. e. dttaw/dlnui > 10~ 10 at uj = 10 2 Hz 
(assuming that the design goals would actually be achieved), 

1 n2 



4al 



11- -\nz s + z; 3 (g s / gi )- 2 



> 10 



21 



(3.6) 



for z s < 10 9 , and either 



logio — < 
9s 



lo Sio 9i ) lo gio 



or 



log 



10 



9i 
9s 



> 



1 



!ogio 9i logio 2. 



/3 < 3/2 



/3 > 3/2 



(3.7) 



(3.8) 



,3 1? 
for z s > 10 9 . 

It may be useful to list approximate forms of the GW spectral distribution, dQcw/dlnw, 
in different regions of parameter space, which we do in Table 1, 





z s < 10 9 


z s > 10 
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Table 1. Approximate forms of (dVlQ\y / d In Cij) in various regions of parameter space, Cc^ — U)\/jZ s , P — 

-log(g s /g 1 )/\ogz s . 

Actually, if one considers also the amplification of electromagnetic (EM) perturbations in 



this scenario ||14j| , one finds an even stronger bound following from the condition Qem < 1) 
i.e. (<7i/<7s) < z 2 s . Such a constraint washes out completely the allowed region of parameter 
space corresponding to the lower row of Table 1. The condition Q < 1 has to be imposed, 
however, for the validity of the linear approach, but it could be evaded in a more general, 
inhomogeneous model of background, in which the evolution of fluctuations is treated non- 
pert urbatively. 

In Figure 1 the allowed region of parameter space, corresponding to the possible detec- 
tion of the spectrum appearing in the upper row of Table 1, and compatible with the various 
bounds on the parameters, is plotted by taking g% = 1 as a reference value. Also shown are 
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the parameter intervals in which other detectors may be useful, for the range corresponding 
to the upper left corner of Table 1. 



Log 10 (G) s /Hz) 




L °gl(Ks 



Microwave cavities Interferometers 
Figure 1. The allowed region in {z s: C/ s ^ parameter space corresponding to the first row in Tclble 1 is the shaded 
region defined by the various constraints. The dashed lines mark the regions in which various detectors may be useful, for the 

range z s < 10 9 and 13 < 3/2. 

We now turn to discuss in more detail the observability of our spectrum as function of 
frequency, limiting our attention, for sake of simplicity, to the frequencies leaving the horizon 
during the dilaton-driven phase [eqs. (|3.1|) , fl3.3|) ]. This spectrum is pictorially described in 
Fig. 2 for the case (3 < 3/2, in terms of the quantity \5h^\ (denoted h c in which 
represents the characteristic amplitude of a stochastic background. The odd-shaped region 
in Fig. 2 shows detection sensitivities for the so called "Advanced LIGO" project, in 
terms of the quantity h 3 / yr defined as the amplitude necessary for detection of a stochastic 
background at the 90% confidence level in a 1/3 of a year (see || for exact definitions). In 
Fig. 2 we can observe clearly that larger amplification goes together with larger red-shift 
for this region of parameter space. For a given red-shift z s , the higher amplitudes, for all 
regions of parameter space, are at higher frequencies. 

In addition to interferometers and bars, microwave cavities may be operated as gravity 
wave detectors for the high frequency range 10 6 — 10 9 Hz. For the MHz range specific 



suggestions [21], [22|] were actually implemented P3| . As can be seen from Fig. 2, the required 
sensitivity for detection of gravity waves in the MHz region is \Sh^\ ~ 10~ 26 , corresponding 
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to h3/ yr of the same order. We leave to experts to study whether or not such a sensitivity 
is accessible with current technologies and may be reached in a near future, with microwave 
cavity detectors or with other experimental apparatus. 




io 2 io 6 io 10 

Frequency in Hertz 



Figure 2. The characteristic spectral amplitude of gravitational waves \5h | . The solid lines show several individual 

II X X-\ TTICLX I 
. The thick dashed line shows the maximum amplitude 1 fl | as a function 

of Z s for Cj s = 1. The dashed lines are lines of fixed g s and therefore lines of constant energy density. VIq\\t is roughly 

the maximal amount of gravitational energy density at a given value of Q s . Also shown in the figure is the odd-shaped region 

marking the sensitivity goals for the detection of a stochastic background according to "Advanced LIGO" project. 

4 Conclusion 

We showed that a rather generic string cosmology scenario leads, naturally, to the production 
of an amplified quasi-thermal spectrum of gravitons during the dilaton-driven phase. This 
spectrum is very stable under modifications of the background and in particular under 0(d, d) 
transformations. The slope of the spectrum may change for modes crossing the horizon in 
the subsequent string phase, but remains in general characterized by an enhanced production 
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of high frequency gravitons, irrespective of the particular value of the spectral index. 

We showed, in particular, that it may be possible to detect such a relic GW background 
with large interferometers for a range of the two parameters characterizing our class of mod- 
els. We would like, however, to encourage the study and the developments of gravitational 
detectors with enhanced sensitivity in the high frequency, KHz - GHz, range. This frequency 
band should be in fact all but a "desert" of relic gravitational radiation that one may expect 
on the grounds of the standard inflationary scenario or from ordinary astrophysical sources. 
Our string cosmology scenario is unique in predicting a strong signal in this range of fre- 
quencies. In general, a sensitivity of Q ~ 10~ 4 — 10~ 5 (which is not out of reach, in the KHz 
region, for coincidence experiments between bars and interferometers @), could be already 
enough to detect a signal, so even a null result at that level of sensitivity would already con- 
strain in a significant way the parameters of the string background, while detection would 
provide a first glimpse at some new and exciting Planckian physics. 

As we stressed, many of our results are independent of details of the string cosmology 
scenario. However, it would be worthwhile pointing out again that, although some ideas 
have been put forward [23J], a solid string-theoretic treatment of the stringy phase, which 
we propose as the necessary bridge between the dilaton-driven and the standard decelerated 
era, does not yet exist. Recent progress ||24| on the interpretation of singularities in string 
theory, as simply a failure to describe physics in terms of the original set of massless fields, 
may shed light on the resolution of this issue. The understanding of singularities in string 
theory would certainly help putting our string cosmology scenario on a firmer basis, and 
may even provide a framework for the calculation of the parameters g s and z s . 
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